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A KINEMATICALLY CONSISTENT TWO-POINT CORRELATION FUNCTION 


J. R. RISTORCELLI* 


Abstract. A simple kinematically consistent expression for the longitudinal two- point correlation func- 
tion related to both the integral length scale and the Taylor microscale is obtained. On the inner scale, in a 
region of width inversely proportional to the turbulent Reynolds number, the function has the appropriate 
curvature at the origin. The expression for two-point correlation is related to the nonlinear cascade rate, or 
dissipation £, a quantity that is carried as part of a typical single-point turbulence closure simulation. Con- 
structing an expression for the two-point correlation whose curvature at the origin is the Taylor microscale 
incorporates one of the fundamental quantities characterizing turbulence, e, into a model for the two-point 
correlation function. The integral of the function also gives, as is required, an outer integral length scale 
of the turbulence independent of viscosity. The proposed expression is obtained by kinematic arguments; 
the intention is to produce a practically applicable expression - in terms of simple elementary functions - 
that allow an analytical evaluation, by asymptotic methods, of diverse functionals relevant to single-point 
turbulence closures. Using the expression devised an example of the asymptotic method by which functionals 
of the two-point correlation can be evaluated is given. 

Key words, two-point correlation functions, turbulence modeling, functionals 

Subject classification. Fluid Mechanics, Aeroacoustics 

1. Introduction. In single-point turbulence closures one is often in the position of needing to approx- 
imate functionals - integrals of functions of the two-point correlation function - Batchelor [1], Kraichnan 
[2], Ristorcelli [3], In aeroacoustical developments in which acoustic radiation due to turbulence is related 
to the statistics of the turbulent source field one is in a similar position: the acoustic radiation is a function 
of diverse two-point integrals, Proudman [4], Ribner [5], Lillcy [6]. 

Batchelor [1], for example, in order to finish his development for the pressure variance in an isotopic 
turbulence needed to evaluate the two functionals 1(f) = f x(f') 2 dx and 1(f) — f x~ 1 (f r ) 2 dx. The function 
f(x) is the two-point longitudinal correlation function. Lilley [6] in an investigation of noise radiated from 
isotropic turbulence required an evaluation of 1(f) — J x 4 (f f ) 2 dx. Proudman [4] in his statistical application 
of Lighthill’s acoustic analogy required approximations to integrals such as 1(f) = f f(f m -\-4x~ l f”— 4r“ 2 /']. 
Ristorcelli [3] in a treatment of weakly compressible turbulence required integrals of the form 1(f) = f xfdx , 
and 1(f) = f xf'(ff f " + 4x l ff rr + Sx 1 f f f f — 4x 2 ff f )dx. These functionals typically involve higher-order 
derivatives of the two-point function. Practical estimates of these functionals are usually made by assuming 
simple exponential or Gaussian behavior for the two-point correlation. It is argued that the functional is an 
integral and therefore only nominally sensitive to form of the assumed two-point function. Unfortunately 
this is not the case when the higher order derivatives appear in the functionals. 

Given the practical need of evaluating these functionals a few attempts have been made to devise, in terms 
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of elementary functions, a useful approximation to the two-point correlation. Several examples of attempts 
are given in Hinze [7], Frenkiel [8], [9], Townsend [10]. Hinze [7] has several amplifications on these issues. 
Attempts have also been made in spectral space: one assumes a spectral energy function shape consistent 
with a Kolmogorov spectrum and transforms to physical space to get the two-point correlation. One obtains 
fractional K n type Bessel functions that are consistent with the Kolmogorov inertial range scalings for 
the structure functions, von Karman [11], and Sirovich, Smith and Yakhot [12]. Such developments are 
focussed at reproducing the intermediate field scalings consistent with the ideal Kolmogorov spectrum: they 
accurately capture the r 2y/3 Kolmogorov inertial range result, Batchelor [15], Monin and Yaglom [16]. 

Our interest, however, is in non-ideal, anisotropic, strained, inhomogeneous, engineering turbulence for 
which an inertial range with an ideal Kolmogorov scaling is unlikely. Such flows are typically computed with 
single-point moment closures in which one parameterizes the turbulence in terms of the Reynolds stresses, 
( v i v j)> an d the dissipation, e. For calculations of such engineering problems, to which such a two parameter 
characterization is applicable, the details of the spectrum are not relevant. Pertinent to this is the fact that 
in physical flows the two-point function is notoriously insensitive to the flow type, Townsend [10], Frenkiel 
[13], [14]. All that appears necessary is that the two-point correlation capture the integral and the dissipation 
length scales. Associated with a characterization of the turbulence by {viVj), and e are two important length 
scales: the integral length scale, £, and the Taylor microscale, A. 

The Taylor microscale is the radius of curvature of the two- point correlation at the origin. The Taylor 
microscale can be related to the viscosity and the dissipation. For a stationary turbulence the dissipation rate 
is the spectral cascade rate. Thus constructing an expression for the two-point correlation whose curvature 
at the origin is the Taylor microscale incorporates one of the fundamental quantities, the cascade rate, 
characterizing turbulence into a model for the two-point correlation function. 

The Taylor microscale is also related to the turbulent Reynolds number, Tennekes and Lumley [19], and 
as consequence the Reynolds number dependence of the curvature at the origin is built into the expression 
for the two-point correlation. Thus the functionals are no longer pure numbers but functions of the turbulent 
Reynolds numbers, / = J(/(x; ^t))» as one might expect. Any satisfactory approximation of the higher- 
order derivatives in the vicinity of the origin, where they are large and a function of the turbulent Reynolds 
number is crucial for any useful evaluation of such functionals. What is desired is a simple analytical 
expression that will be a good approximation for the two-point correlation and its higher derivatives. It 
must also be analytically tractable in terms of elementary functions so that functionals as complicated as 
those given above can be evaluated. Several attempts at such an empirical function are known: Frenkiel 
[8], Townsend [10], von Karman [11], Hinze [7]. None of these attempts satisfy the requirements that allow 
application to situations described in Ristorcelli [3], Ribner [5] Proudman [4], Lilley [6] in which the behavior 
to the two-point correlation near the origin is necessary. 

The purpose of this note is to 1) develop a kinematically appropriate form of a two-point correlation 
function, 2) that incorporates the effect of the turbulent Reynolds number and 3) that is analytically tractable 
in terms of elementary functions. This is accomplished in the following sequence: In §2 the properties of the 
two-point correlation are reviewed. Relationships from statistical fluid mechanics indicating the dependence 
of the function on the Reynolds number and dissipation are summarized. The problem is recognized as a 
two-length scale problem. In as much as there is a large disparity in the two length scales the problem is 
further recognized as one in which there exists a small parameter. Given the existence of a small parameter a 
form of the two-point correlation satisfying all the properties delineated in §2 is proposed and investigated in 
§3. In §3 the boundary layer nature of the problem is used to produce a relationship between the parameters 
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that yields a viscosity independent first integral. The details of the calculations are provided as an example 
of the asymptotic method of evaluating similar functionals. 

2. Mathematical and observed properties. The diverse kinematic and observed properties of the 
two-point correlation are now given. The relationship between curvature at the origin, Taylor microscale, 
dissipation, and Reynolds number is summarized. The definition for the correlation function is obtained from 
< Vi(x)Vj(x + r ) >— < ViVj >= ^ Rij(r ), where k = \ < VjVj >. The isotropic portion of the two-point 
correlation can be written as 

(1) %(r) = -^/'+ (/+lr/')f« 
where f = f(r) the longitudinal correlation is defined: 

2k 2k 

(2) < vi(0)vi(r) >=< vivi > f(r) = y/(r) = yi?n, 

von Karman and Howarth [17], Batchelor [18], Monin and Yaglom [16]. Attention will be restricted to the 
isotropic component of the longitudinal two-point correlation correlation. The correlation function, /(r), 
has the following properties 

m = 1 

/'( o) = o 
r(o)^-i/A 2 

/(oo) -+ 0. 

The function is also even: / 2n+1 (0) = 0. It is typical to work in nondimensional variables: one rescales 
according to x = rji where £ is the integral length scale of the turbulence. The integral length scale is 
defined as | ki= J 0 °° < v\v[ > dr and thus 

rOC 

(3) / f(x) dx = 1. 

Jo 

The curvature at the origin, in nondimensional coordinates, is then specified as f"(x = 0) = —£ 2 /\ 2 . The 
Taylor microscale, A, can be shown to be equal to the spectrally weighted length scale 

J_ _ JjfE{K)dK 

' ’ A 2 jE{n)dK ' 

The Taylor microscale is an intermediate length scale; smaller then the outer integral length scale and larger 
than the inner Kolmogorov scale, Tennekes and Lumley [19]. 

The turbulent Reynolds number is defined as R t = u c ijv where u c is a characteristic fluctuating velocity, 
u c = (I&) 1 / 2 . The Kolmogorov scaling e = akU^/£, [19], is used to eliminate the length scale: thus R t = 
ak^k 2 /{ve). Note that this is the traditional Reynolds number, Tennekes and Lumley [19]; the Reynolds 
number definition used in contemporary DNS is a factor 9 larger and does not realistically reflect the relative 
magnitude of inertial to viscous forces. The isotropic portion of the dissipation tensor, €ij = v < Uj,k >, 
can be written in terms of the two-point correlation tensor: 

2 2 2 

(5) €jj = 2e == v < Ui,k Uj ik >— -kvRjj^kk |o = (0) = — fci/15/A , 

Batchelor [15]. The curvature, A, has been related to the spectral cascade rate of energy, e : 

(6) e = bkv/X 2 , 
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Tennekes and Lumley [19]. The cascade rate, e, is a measure of one of the most fundamental characteristics 
of turbulence its nonlinear decorrelating effect; any model for a two-point correlation must be dependent 
on this quantity. The curvature, A, is also related to the Reynolds number. The Reynolds number, Rt = 
^k 2 /(i/e), with otk — 1 can be used to eliminate e. One finds, Tennekes and Lumley [19], that 

(7) t 2 /X 2 = ±R t = f"(x = 0). 

The dependence of the curvature at origin to turbulent Reynolds number needs to be embedded in any 
empirical expression of the correlation function. Moreover any model for the two-point correlation, / = 
/(x; Rt), must have 

rOC 

(8) / f(x; R t ) dx= 1, 

Jo 

a first integral independent of Reynolds number (or, equivalently, viscosity). 

An additional constraint for any model of / can be found from the problem of the final period of the 
decay. In the final period of the decay, first treated by von Karman and Howarth [17], (see also [18], [7]), 
the two-point correlation is found to be a Gaussian: as R t — » 0 

r 2 

( 9 ) / ~ e 9t ' t , 

Note that in the linear decay problem &vt = A 2 ~ l/Rt and 8ut — > oo corresponds to Rt — > 0. This function 
has the appropriate properties at the origin: it is an even (differentiable) function with finite curvature 
specified by the Reynolds number. The Gaussian, however, cannot satisfy J 0 °° f(x ) dx = 1 in a way that is 
independent of Reynolds number. 

It is empirically observed that the decay of the correlation function in high Reynolds number turbulence 
can be usefully approximated as / x , Hinze [7], Frenkiel [8], Townsend [10]. A similar behavior for the 
two-time correlation is discussed in Pope [20]. In as much as e~^ is an easily integrable function it is used 
in estimates of integrals of the two-point correlation. It does not have any of the required properties near the 
origin. This, nonetheless, has not stopped its application as an estimate in a number of situations, [2], [7], 

(10) . Its application being understood to be limited to issues related to larger scales of the motion, [20]. The 
form cannot, of course, be used in any functionals that feature a dependence on higher order derivatives 
which are very large near the origin and cannot be represented by e - ^, which is non- differentiable at the 
origin. 

3. Mathematical representation. Many of the empirical functions chosen to approximate the two- 

point correlation in the estimation of functionals of / are single parameter curves like e~°^ or e~~ ax which 

while satisfying one constraint do not reflect the two-length scale nature the correlation function. For 

example, near the origin, / ~ e ~ ^ has positive curvature while the real correlation function has negative 

curvature. The use of / ~ e ax which has proper sign for the curvature will not have the proper magnitude 

of curvature at the origin and produce, as is required, a viscosity or Reynolds number independent integral 

2 

length scale. Furthermore / ~ e _ax is inconsistent with the experimentally the observed exponential decay, 
e -a|x|, jgj Yet the choice / ~ e~^ while a useful approximation to the experimental data, Ftenkiel [8], 
Townsend [10], has finite (discontinuous) slope at the origin and is an odd function. In addition the positive 
curvature of / ~ e~ ax at the origin which is unacceptably inconsistent with the dissipative nature of the 
small scales. Any single parameter two-point correlation model cannot capture both of the length scale 
properties of a turbulence correlation: the integral length scale, and the Taylor microscale, A. 
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A simple two parameter expression for /(x; Rt) in terms of elementary functions satisfying the A and £ 
properties is now considered. The properties that any model for the two-point correlation must satisfy are: 

• The two-point function has the properties: /( 0) = l,/'(0) — 0,/(oo) — > 0 and f 2n+l (0) = 0. 

• The curvature at the origin is specified /"( 0) = — 1/A 2 . 

• Its outer length scale is specified by the normalization / 0 °° f(x)dx = 1. 

r 2 

• It satisfies the exact result for the low Reynolds number limit, / ^ e~ ** . 

• It has, as is observed for large Reynolds number, an exponential decay for large x, / ~ e~ x . 

We are now in the position of proposing a function that has all these properties: 

(10) f{x\€,b) = 


The function has two scale parameters, e and 6, that will be used to satisfy the two length scales constraints. 
The two asymptotic forms of /, in as much as they highlight the two scale parameters, are worth considering. 
For small x 


(ii) 


0, 


_ii i 

> c (l + b 2 x 2 /e 2 )) 1 /2 



e c 


For a high Reynolds number turbulence, the quantity e, will be seen to be a small parameter: c ~ 1 /Rt 
number. The small parameter, e, forms a boundary layer in the vicinity of the origin, x << e 2 , and the 
function has the required large and positive curvature. For low Reynolds number, e ~ 1 and the Gaussian 
behavior of the von Karman and Howarth result is obtained. For large x the function is written as 

(12) x — ► oc, / ~ e ^ (i + e 2 /(& 2 * 2 )) 1 / 2 e~%. 


The values of the two scale parameters are now related to quantities describing the turbulence. 

It was shown above that the curvature of the correlation function at the origin must satisfies /"( 0) = 
— ~R t . Differentiating / twice, /"( 0) = — 2 = — j$Rt and the small parameter is determined by the 
turbulent Reynolds number: 


(13) 


€ — 


15 

TV 


The two-point correlation also satisfies the normalization condition, 


(14) 


h 


(/(«, b) = r 

Jo 


/(x; e, b) dx — 1 


which will serve to specify b in terms of Rt in a way that I\ is independent of Rt . The integral I\ is not 
tractable but an asymptotic analysis produces simple analytical results. 

The integral I\ is an example of integral that has local and global contributions. The usual Laplace 
method of evaluating such integrals with large parameters is not applicable, however, the boundary layer of 
the expression for /(x; R t ) can be exploited. Near origin for x < e there is a boundary layer whose width 
scales with 1 fR t . In the boundary layer region / does not change appreciably and the local contribution 
to the integral I\ scales with the width of the region, of O(e). The global contribution to /i occurs over a 
region 0(b) for x > e. In this region, / ~ e~ x and / 0 °° e~~ x dx = 1 and the global contribution is 0(1). With 
these ideas in mind the the interval of integration is subdivided into local and global regions in which the 
two-point function has, respectively, the near and far field behavior given by (11) and (12): 

nS rOO 

(15) I u + /loo = / / (x) dx + / f{x)dx — 

Jo Js 
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Here 5 satisfies 6<<J<<6<<00 but is otherwise arbitrary. The fact of its arbitrariness can be used to 
validate the success of the asymptotic evaluation of the two integrals I\ t and the combination I u + iioo 
must be independent of 5. As the contribution to each integral is proportional to the their intervals one 
might expect expect 7i c ~ 0{8) and 1^ ~ 0(6); since ^ << 1 the major contribution to I\ is from I\ ^ 
and one might expect 6 ~ 1 to leading order. The local and global contributions to the integral I\ are now 
evaluated. 


Inner integral. For the inner region the change of variable x — er} is made. A Taylor series expansion 
for the exponential, since e < < 1 , is used to produce the following expression valid for x 0(e), 


(16) 


f Q =e (l + & 2 rj 2 )) 1/2 % 1 


where 0 < rj < 8/t. The integral I\ t is rewritten 


Uc 


fS/e 

L 1 


1 


(1 + 6V) 1/2 


(1 + 6 V )) 1/2 


+ 0(e 2 ) ] dr) 


+ 0(e 2 ) 


(17) 


le6S 2 b 2 1/2 

= S ~2V^ + 1) 


1 e 2 bS 

+ 5P in l7 + ( -^ + 1 )/ l + • 


Re-expressing the terms as a function of ^ and realizing that both e << 1 and | << 1 and, using the Taylor 
series expansions for the quadratic and the logarithm, one obtains 


(18) 


_ r 1 6 2 . 1 e 2 , If 2 , r 26<5, 4 , 

Iu ~ 5 ~2J {1 + 2 + 2P M V ] + ° {e 


Outer integral. For large x the two-point function is rewritten 
(19) /= e _ ^ {i+. a /c62 3 ,a ) )i / 2 _ e -f [i-ijif? + |(^t ) 2 + -I _ e -f 


1^; 
e 2 J 5 


_ 12. ( _ 
e 86 ' 6 


where we have used the binomial expansion since x > 5 and therefore << 1. The last two factors are 
expanded using the Taylor series for the exponential. To 0(e 2 ) one obtains 

rOC -I 9 7*00 — ^ 

(2°) I X co = l e-*dz + j s —dx. 

The first integral is I\oo = be~i and is to be expanded in powers of | since | << 1. The second integral is 
the exponential integral. The exponential integral is also expanded in powers of | : 


( 21 ) 


le 2 f°° e~% , 

tPJ. — dX 


If f_ 

2 b 2 7 


, ,6. 6 IS 2 1 S 3 

ln ty + b ~ lb 2 + 18 b 3 ' 


where 7 « 0.57721 is Euler’s constant. The integral becomes 


( 22 ) 


1.6 


le 2 


/roc = 6[1 - r + r(r) + - 1 + 53 ["7 


, A S 

Zn( 6 )+ 6 


b ' 2 v f> / ' ' J ' 2 i> 2 

Using the expressions (18) and (22) in I\ f + J\ ^ = 1 produces, to leading order, 


If J_f 

4 6 2 18 b 3 


(23) 


1 c 2 1 b 2 

6 = 1 + 2ph-2i-'"< 2 7>I- 


Note that the expression does not depend on the arbitrary scale factor J; this is a vindication of the procedure 
and can be used to check for errors. The nonlinear expression for b is solved iteratively. Only one iteration 
is required. Setting b = 1 in the right hand side produces 


(24) 


6=1 + ^ 2 t 7 -l"^)l- 
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To a very good approximation b ~ 1 for high Reynolds (as low as Rt > 100). The expression for the two-point 
correlation, with b m 1, becomes 

(25) /(x; R t ) = e _ «?i 5 +«P) T7 \ 

This, (25), is the expression for a kinematically consistent two-point correlation whose curvature at the 
origin is given by the Taylor microscale. The small and large x limits are, respectively / — + e ~rs RtX and 
f — > e ~ x . The asymptotic evaluation of ii was used to obtain b = b(Rt) such that I\ was independent 
of Ru b » 1 was indicated to be a good approximation. The asymptotic evaluation of the I\ integral 
was described in some detail: it is an illustration of similar asymptotic procedures needed to evaluate the 
functionals described in Batchelor [1], Proudman [4], Ribner [5], Lilley [6], Ristorcelli [3]. 

The fact that b w 1 is a result of the fact that I u ~ 0(e) « hoo- This will only be the case for 
functionals that do not include higher derivatives. Higher order derivatives and their products will scales 
with R™, n > 0, near the origin and the fact that the interval scales with R^ 1 will mean that the local 
contribution will be larger than the global contribution, 7 C >> Toe- 

4. Discussion. The kinematic two-point correlation has been obtained with the idea of applying it to 
functionals calculated for non- ideal flow situations in which a Kolmogorov spectrum is not expected. To 
this end we have not incorporated the r 2 / 3 scaling, [15], associated with an ideal Kolmogorov turbulence. 
It is, at the expense of convenience and simplicity, possible to incorporate the Kolmogorov behavior. [It 
involves another free parameter and a length scale.] Due to the complicated expression needed to build 
in Kolmogorov behavior, the fact that information necessary to specify the additional length scale can not 
be obtained from turbulence closures, and the nonideal nature of practical flows, consistency with the r 2 / 3 
scaling has been foregone. 

The result, (25), is more general than if we required consistency with the Kolmogorov behavior. It 
becomes an interesting kinematically consistent expression for any random velocity field with distinguished 
inner and outer [integral and dissipation] scales. It is straightforward to produce similar expressions for 
either Eulerian or Lagrangian two-time correlations, [7]. In a Langevin approaches to turbulence, Pope [20], 
two-time correlations of the form p(s) = e~ 8 ^ r are used; higher order derivatives are again, unphysically 
discontinuous at the origin. The model, e -S//T , is the correlation function of the Orhnstein-Uhlenbeck process 
(colored noise) with non- differentiable velocity increments for which Pope [20] has a nice discussion relevant 
to the turbulence problem. Should statistics to which the differentiability property or the low Reynolds 
number limit be important than an expression of the form (25) is straightforward. The boundary layer at 
the origin is a temporal layer and the same procedure with two relevant time scales, the Lagrangian integral 
time scale and the Lagrangian dissipation time scale (Hinze [7]), is applicable. 

5. Closure. A simple kinematically consistent two-parameter expression for the two-point spatial cor- 
relation function has been obtained. The model two-point expression has curvature at the origin given by 
the Taylor microscale and satisfies the integral constraint associated with a Reynolds number independent 
integral length scale. In as much as the curvature at the origin is related to the dissipation the model 
constructed incorporates one of the fundamental quantities characterizing a stationary turbulence, e. As 
the curvature at the origin can also be expressed in terms of the turbulent Reynolds number, a measure of 
nonlinearity, the model expression for the two-point correlation includes the decorrelating effects associated 
with the [nonlinear] cascade. 

The expression is meant to be used to provide estimates, using an asymptotic procedure described, 
of functionals used in single-point moment models. The expression is purely kinematic; no appeal to the 
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dynamical two-point equations has been made. Any development based on the dynamical two-point equations 
for a general anisotropic inhomogeneous turbulent flow undergoing strain is prohibitively complicated and 
unclosed. In as much as the shape of the two-point correlation is insensitive to flow situation, Townsend 
[10], a kinematically consistent expression will capture the two- length scale features seen isotropic Frenkiel 
[8], strained Rogers [21], and boundary layer Frenkiel and Klebanoff [13], The expression is consistent with 
several kinematic properties including differentiability at the origin and allows one to express two-point 
behavior in terms of two quantities carried in a typical moment closure scheme, the kinetic energy, k } and 
dissipation, e, of turbulence. 
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